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Abstract 


Metallorganic molecules have been proposed as excellent spin filters in molecular spintron- 
ics because of the large spin-polarization of their electronic structure. However, most of the 
studies involving spin transport, have disregarded fundamental aspects such as the magnetic 
anisotropy of the molecule and the excitation of spin-flip processes during electron transport. 
Here, we study a molecule containing a Co and an Fe atoms stacked between three cyclopen- 
tadienyl rings that presents a large magnetic anisotropy and a S=l. These figures are superior 
to other molecules with the same transition metal, and improves the spin-filtering capacities 
of the molecule. Non-equilibrium Green’s functions calculations based on density functional 
theory predict excellent spin-filtering properties both in tunnel and contact transport regimes. 
However, exciting the first magnetic state drastically reduces the current’s spin polarization. 
Furthermore, a difference of temperature between electrodes leads to strong thermoelectric ef¬ 
fects that also suppress spin polarization. Our study shows that in-principle good molecular 
candidates for spintronics need to be confronted with inelastic and thermoelectric effects. 


Introduction 

Molecular spintronics is a thriving field driven by advances in shrinking electronic devices using 
molecules^ and by the extraordinary properties of spin transport . 2,3 Not only are molecules com¬ 
plex enough to attain dedicated functionalities, but they are identically replicated and cheap to man¬ 
ufacture using chemical synthesis. Recently, it has been possible to address individual molecules 
while taking advantage of their hierarchical growth to create structures of increasing complexity. 4 - 
Molecules can become fundamental pieces of the ever shrinking device technology. 5 Additionally, 
molecules show a great diversity of magnetic properties that can be successfully tailored, such as 
spin-crossover molecules,- molecular magnets,- spin-filtering molecules , 2 molecular spin valves- 
and molecular switches . 7 Molecular spintronics is then a rich field which promises scientific and 
technological breakthroughs. 

An interesting functionality that has been sought after in molecules is the capability of select- 
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ing one spin to be transmitted in a given spintronic device. 9 In order to achieve this, the molecule 
presents spin-polarized frontier orbitals, with one of the spins more coupled to the contacting 
electrodes. In this way, the coupled molecular orbital has a larger contribution to electronic trans¬ 
port, favoring the transmission of one spin species. Typically the molecular spin polarization is 
achieved by using complexes where the metallic atom (or atoms) present an open-shell configura¬ 
tion. The ligand field of the rest of the molecule on the metallic atoms lead to interesting physics: 
different spins can be present within the small energy scale of the ligand field. 6 This is particu¬ 
larly true in the case of a sizable spin-orbit coupling (SOC) and a spin larger than 1/2, because 
the ligand field creates a magnetic anisotropy due to the SOC that can fix the orientation of the 
molecular spin leading to the appearance of molecular magnets. 10 However, even in the absence 
of a fixed magnetic-moment orientation, many different molecular systems have been signaled as 
spin filters, because transport is basically dominated by one of the electron’s spins. Indeed, recent 
works—show that in the absence of a magnetic center, radical molecules can be used leading to 
spin-polarized electron transport. 

Large spin polarizations have been predicted for the family of molecules made from interca¬ 
lated sequences of organic rings and transition metals. Examples of these molecules are benzene- 
vanadium ensembles, 13 benzene-cobalt, 9 cobaltoccne- 4 and ferrocene and 1-D ferrocene-based 
wires. 15 Complete studies of different stacking of cyclopentadienyl (Cp) and transition metals 
(TM) or benzene and transition metals have also been performed. 1617 Stacking two different TM’s 
has been less common. Some calculations suggest that infinite sequences of stacked TM-Cp 
present exotic electronic structure with different magnetic ordering depending on the used TM 
atom. 18 Here, we propose a new molecular spin filter by stacking an iron and a cobalt atom be¬ 
tween three cyclopentadienyls (Cp-Fe-Cp-Co-Cp).We performed non-equilibrium Green’s func¬ 
tions (NEGF) simulations to evaluate the transport properties of this molecule, CoFeCp 3 , based 
on density functional theory (DFT). As expected, the spin polarization obtained in transport ap¬ 
proaches 100%. Moreover, the hybrid magnetic structure of this molecule leads to a ferromag¬ 
netic coupling between the magnetic centers, where most of the magnetization is localized on the 
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cobalt atom. Due to the sizable spin-orbit coupling of Co, the Cp ring induce a sizable magnetic 
anisotropy energy (MAE) which is very interesting for spin-filtering applications. However, this 
same energy scale sets the energy scale for the first spin excitations that can drastically reduce the 
spin-polarization in the electron current. 19 20 We evaluate here the effect of bias in reducing the 
spin polarization as spin-flip processes become energetically accessible. 

Our calculations show that transport takes place through the molecular electronic structure 
based on its ^-orbitals. The broken-symmetry electronic structure of CoFeCp 3 leads to frontier 
orbitals of different nature and spin. In contact with metallic electrodes, only the tails of the reso¬ 
nances caused by the molecule-electrode interaction contribute to transport. Hence, transmission 
changes rapidly with energy near the Fermi energy which should lead to large thermoelectric ef¬ 
fects. 2 ^ 23 Moreover, the thermoelectric properties should be different per spin, which can lead to 
spin currents even in the absence of charge currents.—“ 2 - 

Theoretical methods 

In order to perform the calculations of this work, we have mainly used two density-functional the¬ 
ory (DFT) packages. VASP— - — has been used to explore the adsorption of the CoFeCp 3 molecule 
on the Cu(l 11) surface and also its magnetic anisotropy. Geometrical effects when a second elec¬ 
trode (another Cu(lll) surface) was approached, have been evaluated with VASP. However, the 
bulk of the calculations has been performed using the Siesta package. 33 These calculations con¬ 
firmed the results obtained from VASP and permitted us to perform electronic transport calcula¬ 
tions using TranSiesta.— 

We optimized the structure of the CoFeCp 3 lCu(l 11) interface, using density functional theory 
(DFT) at the spin-polarized generalized gradient approximation (GGA-PBE) level, as implemented 
in VASP. 27-32 In order to introduce long-range dispersion corrections, we employed the so called 
DFT-D2 approach proposed by Grimme. 35 We used a plane wave basis set and the projected aug¬ 
mented wave (PAW) method with an energy cut-off of 400 eV. A 19-A thick vacuum region was 
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used to decouple the surfaces of consecutive slabs in the supercell approach used in VASP. The 
surfaces were modeled using a slab geometry with five Cu layers and a 3 x 2y/3 unit cell. Such an 
unusual unit cell have been chosen based on experimental data for ferrocene (FeCp 2 ) molecules, 
which can be seen as one of the building blocks for CoFeCp 3 . Self-assembled monolayers of 
ferrocene shows a 6 x 2y/3 periodicity with two molecules per unit cell. 36 Published calculations 
yield that these two molecules do not interact between them. 37 Therefore, we decided to carry out 
our calculations using a smaller 3 x 2\/3, which is still large enough to prevent interactions among 
adsorbed molecules. 

During the geometry optimizations, we allowed for the relaxation of all atoms of the molecule 
and of the two-topmost layers of the Cu surface until the atomic forces were smaller than 0.02 eV/A. 
A 7x7x1 k-point sampling of the first Brillouin zone was performed using the Monkhorst-Pack 
method. 38 

Transport calculations were carried out from first-principles with a method based on nonequi¬ 
librium Green’s functions (NEGF) combined with DFT as implemented in the TranSiesta pack¬ 
age. 34 The open-boundary system is divided in three distinct regions breaking the periodicity along 
the transport direction. The central part is the scattering region and the other two regions are the 
semi-infinite left and right electrodes, formed by periodically repeating six layers of bulk copper. 

The most favorable configuration after geometrical optimization of the CoFeCp 3 lCu(l 11) inter¬ 
face was used to build the scattering region. As illustrated in Fig.Q] the scattering region was com¬ 
posed of one CoFeCp 3 molecule connected to two Cu(l 11) surfaces, left and right, each formed by 
8 active layers of a 3 x 2y/3 cell. It is important to stress that two-probe system geometries were 
obtained after geometry optimizations using VASP. Dispersion corrections were described through 
the semi-empirical DFT+D2 scheme. Hence, the role of dispersion forces on the transport results 
is implicitly considered trough the optimization of the junction’s geometry. 

For transport calculations, the valence electrons wave functions were expanded in a basis set 
of local orbitals. A double-£ plus polarization (DZP) basis set was used to describe the molecu¬ 
lar states and and a single-£ plus polarization orbitals (SZP) basis set for the copper electrodes. 
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Diffuse functions were also included to describe surface electrons. The use of a DZP basis set 


to describe the molecular states is mandatory in order to yield correct transmission functions. In¬ 
deed, a SZP basis set led to a shift of the main molecular peaks of ~ 0.3 eV with respect to the 
DZP molecular peaks. However, using a DZP for the full system does not alter the transmission 
functions noticeably. Therefore, the chosen basis set seems to be a good compromise between 
computational cost and quality. We employed the GGA/PBE functional— and norm-conserving 
Troullier-Martins pseudopotentials. 40 A 11x11 in-plane k-point mesh was adequate to obtain suf¬ 
ficiently accurate transport results. 

The spin-polarized electron current I G (<7 =j",|, denoting majority a minority spin channels 
respectively) was calculated using the Landauer-Buttiker expression:— 



( 1 ) 


where T a (e, V) is the transmission function for an electron of energy £ and spin <7 when the bias 
voltage between the two electrodes is V. In eq.??, f(e,ji v ,T v ) — (1 + exp(£ — ii v )/k B T v )~ l is 


the Fermi Dirac distribution of electrode v (v = L,R, left and right electrodes respectively) with 


temperature T v and chemical potential /i v (note that V = {jl F — /i F )/e). The electron charge is 
given by e and Planck’s constant by h. 

In the linear-response regime, I G can be approximated as— 


I a ~G a V + G a S a (T L -T R ) 


( 2 ) 


where G a and S G are the spin-dependent conductance and Seebeck coefficient which are calculated 
at zero bias voltage (V = 0) as 


Ga = —Kq g (Efi T), 
n 


(3) 


and 


1 K lG {E F J ) 
e\TK 0G (E Fl Ty 


(4) 
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where Ep = jl; = jlp is the Fermi level and 


K a . a (E F ,T) = -J df ( £ f e F ’ I " } ( £-E F ) a T ff (e,0) de 

with a = 0,1. The total electronic conductance is given by G = G| + G|. 

Finally, the spin-filtering capabilities of the molecular junction is analyzed in terms of the spin 
polarization of the current, CP, defined as 


CP=(/ t -4)/(/ t +4)xlOO. (5) 

When both the temperature difference and bias voltage between left and right electrode are zero 
(i.e. V = 0 and T l — Tr — 0) the spin-filtering capacities are evaluated using the spin polarization 
of the transmission function at the Fermi energy. The corresponding quantity is called spin-filter 
efficiency 42 ' 43 and is defined as 


SFE = (T t (E Fenni , 0) - T \{E FermU 0))/(t ^{Epermi, 0) + T \{Ep ermU 0)) X 100. (6) 


Results and Discussions 

In this section, we analyze and discuss the results obtained for CoFeCp 3 as a spin filter in the 
transport of electrons between two copper electrodes. The section is divided in several subsec¬ 
tions to give a thorough view of the properties of this molecular device. The first subsection 
analyzes the isolated molecule and compares it to related molecules, explaining why CoFeCp 3 is a 
good candidate for a spin-filter device. The second subsection analyzes the adsorbed molecule on 
Cu(l 11). The third subsection is devoted to electron transport in the elastic regime in the absence 
of thermoelectric effects, both for tunneling and high-conductance regimes. The modification of 
the spin-filtering capacities when spin-flip processes are allowed is evaluated in the following sub¬ 
section. This section is finished by a detailed account of the effect of thermoelectric effects in the 
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properties of CoFeCp 3 as a spin filter. 


Gas-Phase CoFeCp3 

As shown in Fig.[T] we considered two types of initial structures for CoFeCp 3 molecules: eclipsed 
and staggered (D 5 /, and symmetries, respectively). In agreement with previous results obtained 

for ferrocene, FeCp 2 ,“ the eclipsed conformer is slightly more stable than the staggered one (the 
computed energy difference is 58 meV). 

In both conformers, the ligand field splits the degenerated Co/Fe (TM) d levels into one d v 2 
(a\) and two doubly-degenerated d xy — d x i_ y i (ei) and d zx — d yz {e\) orbitals. Depending on their 
symmetry and energy position, these orbitals mix to a different degree with 2p states of the C 
atoms. For instance, the highest occupied molecular orbitals (HOMO), for majority (HOMOjO 
and minority (HOMO),) spin channels, schematically shown in Fig. [2l have ~ 50% TM-<?| and 
~ 90% TM-t '2 character, respectively. On the other hand, lowest unoccupied orbitals (LUMO) 
for majority (LUMOjO and minority spin channels (LUMOj,) (see Fig. [ 2 ]) present ~ 50% and 75 
% TM-g| character, respectively. This picture agrees well with the ligand-field splitting of the 
c/-clcctron manifold in Ds-symmctry. 

The Cp ligands roughly contain one electron. Hence, the TM atoms approximately are in d 6 
(Fe) and d 1 (Co) configurations, see Table [0 The lowest-energy conformation corresponds to the 
low-spin one, hence filling the ligand-splitted d levels for Fe and Co leads to a spin 1 molecule. 
This is confirmed by our calculations, regardless of the used exchange-and-correlation functional. 
From this picture, we see that Co will host the spin one, and Fe will have spin zero. This is 
in agreement with the zero spin of ferrocene. However, cobaltocene (C 0 CP 2 ) is spin 1/2. The 
difference stems from the presence of a Cp between Fe and Co in CoFeCp 3 . Indeed, CoFeCp 3 
is not a ferrocene plus a cobaltocene. Plotting the spin distribution for CoFeCp 3 , we confirm the 
above results: spin is largely localized on the Co atom, and the Fe atom is basically not magnetic. 

The large spin-orbit coupling of Co, leads to a sizable MAE induced by the Cp’s ligand field. 
We have evaluated the MAE and we obtain that the Co-Fe axis is a hard axis. This means that the 
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magnetic moment of the molecule lies in a plane parallel to the Cp’s. The transversal anisotropy 
is negligible. Hence the magnetic moment is not fixed in a particular direction in the Cp’s plane. 
The MAE is 1.64 meV for both conformers. This is the energy needed to change the magnetic 
moment from the easy plane to the hard axis. Since the magnetic moment corresponds to 5 — 1, 
the molecular ground state is doubly degenerate and formed by the spin components |5 Z | = 1. The 
first excited state is S z — 0. Hence, the magnetic moment will be localized in the Cp’s plane as 
long as the bias between electrodes is not large enough to flip the spin from \S Z \ — 1 to S z — 0 as 
will be discussed below. These results have been obtained in the gas phase and are, in principle, 
not valid for the adsorbed molecule. As we will see in the next section, the molecule is basically 
physisorbed on Cu(lll) without charge transfer or any interaction from the substrate other than 
dispersion forces. Hence, we expect that the gas-phase MAE be a good approximation to the MAE 
of the spin-filter device. 

These data indicate that CoFeCp 3 is a small molecule with an important spin that is fixed 
to a plane contained by the Cp ligands, with a pinning energy (MAE) of 1.64 meV. Hence, the 
molecule can in principle polarize an electronic current to a direction perpendicular to the axis 
of the molecule. It is interesting to compare this molecule with similar molecules. Co2Cp3 or 
Fe 2 Cp 3 will not be good spin filters. 20 The presence of an odd number of Cp leads these molecules 
to present a low-spin configuration 5 — 1/2, which is not subjected to any magnetic anisotropy 
and cannot be molecular magnets. Molecules with an odd number of Cp’s and only one type 
of TM atom such as Fe or Co, will probably not be good spin filters either because they show 
antiferromagnetic coupling with its corresponding 5 = 0 ground state. Infinite chains of CoCp 18 
also show antiferromagnetic ordering and hence a 5 = 0 ground state. The case of FeCp chains is 
more complex. For short molecules, the ground state is the low-spin configuration 5 = 0, however 
as the chain grows larger, a half-metallic ferromagnet develops that can eventually be an excellent 
spin filter.— Here, we propose something simpler, just a CoFeCp 3 molecule. 
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Table 1: Charge distribution and magnetization for the isolated molecule (evaluated with 
Siesta and Mulliken-charge analysis). Total magnetization is 2/ig (S = 1). 


Element 

Total Charge (Mulliken) 

Magnetization (/i#) (Mulliken) 

Fe 

6.687 (d states= 6.191) 

0.513 (d states= 0.475) 

Co 

7.881 (d states= 7.314) 

1.731 (d states= 1.684) 

C 

62.772 

-0.262 (-0.228 Cp in between) 

H 

14.667 

0.018 



(b) 


Figure 1: (a) Schematic top view of D 5d and D 5h conformers within the 3 x 2-\/3 unit cell. 

Yellow: Fe atom, violet: Co atom, orange: C atoms, grey: H atoms. The black lines represent the 
surface unit cell, (b) Lateral view of the scattering region used in transport calculations. 

Adsorption of CoFeCp3 on Cu(lll) 

As a first step we carried out full geometry optimizations for a single FeCoCp 3 molecule with the 
Fe-Co axis initially located on the high-symmetry sites of Cu(lll): top, bridge, hollow-hep, and 
hollow-fee. The most (least) stable final configuration corresponds to the molecule adsorbed on 
the hollow (top) site at an average distance of 2.65 A (2.78 A) from the surface. However, the 
energy difference between top and hollow adsorption sites is only 65 meV. Since the computed 
equilibrium points are spatially very close we do not expect to have large energy barriers between 
the points, leading to an overall small diffusion barrier. 
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Charge population calculations using the Bader scheme^- point to negligible charge transfer 
between molecule and surface. In addition, neither the geometrical structure nor the electronic 
characteristic of the molecule seem to be strongly affected by the adsorption process. As a result, 
the adsorbed molecule maintains its gas-phase electronic and magnetic properties. This is further 
corroborated by a deep analysis of the contributions to the total adsorption energy. 

The main contribution to the adsorption energy, E ads , comes from dispersion forces (E vdw ). 
Indeed, the evaluated adsorption energy on the hollow site is E ads ~ —1.19 eV. The contribution to 
this adsorption energy is mainly due to the van der Waals component, E vd \y ~ —1.28 eV that is re¬ 
duced to the final E ads value by the repulsion with the electronic cloud of the surface. Interestingly, 
if vdW interactions are turned off in the calculations, the molecule feels the repulsive forces and 
reaches an adsorption distance of 3.22 A with a very small adsorption energy (E ads ~ —0.13 eY) 
that is probably not meaningful. Nevertheless, these results show that the molecule binds solely by 
the action of van der Waals forces. 

Transport properties of CoFeCp3 on Cu(lll) 

In the present section, we show the results of our electron transport calculations for a CoFeCp 3 
molecule between two Cu(lll) electrodes with special emphasis on spin filtering. The first re¬ 
sults correspond to the electron transmission across the molecular junction at zero bias. First, the 
tunneling regime is analyzed, where the right electrode is kept at a distance much larger than the 
adsorption one. Then, we analyze the contact regime, also at zero bias. The third subsection ex¬ 
plores bias effect in the more interesting case of the contacted junction. And finally, motivated 
by the slopes of the transmission function at the Fermi energy, we compute the behavior of the 
molecular junction with respect to a temperature gradient and the related thermoelectric effects. 

All the results of this section have been evaluated for the eclipsed (D 5 / ; ) molecular conformer. 
The very similar data about the staggered conformer can be found in the Supporting Information. 
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Transport in the tunneling regime 


Figure [3] (a) shows the transmission spectra at zero bias for a left-(right-)electrode-molecule dis¬ 
tance of d\ =2.65 A (di =5.15 A). We approximate the zero-bias conductance by the transmission 
at the Fermi level, Ep. Hence, the conductance is G(Ep ) = 4.53 x 10~ 3 Go where Go = 2 e 2 /h is 
the quantum of conductance. The small value of G(Ep ) shows that this setup corresponds to the 
tunneling regime. 

Concerning the effect of the electrodes on the geometry and electronic properties of the molecule, 
we conclude that it seems to be very small. Within this geometry i) the structural parameters of 
the molecule are very similar to the ones obtained in gas phase ii) the total charge transfer to the 
molecule is very modest (0.098 e) and iii) the total magnetization of the molecule is 2 jl B where the 
partial contributions coming from Fe, Co, C and H (i.e 0.684 /J, B , 1.567 /i 5 , -0.271 /i B and 0.022 
f.l B , respectively) are close to the values reported in Table I for the isolated molecule 

Figure [3] shows the extraordinary spin-polarization induced by the molecule. The majority-spin 
channel transmission (t \{Ep erm i)) is roughly two orders of magnitude higher than the minority- 
spin one (t \{Ep erm i)). As a result, the SEE given by eq. ?? approaches 100 % (more precisely, 
SFE = 98%). 

The transmission of Fig. [3] (a) implies that transport is mainly determined by the hybridization 
of surface electronic states with the frontier molecular orbitals. To get a deeper understanding 
of the different features observed in the transmission function, we plot the density of states pro¬ 
jected (PDOS) onto the frontier orbitals that we analyzed above, namely, the doubly-degenerated 
HOMO’s and LUMO’s. These PDOS are depicted in Fig. [3] (b). The PDOS peaks perfectly 
match the transmission ones, permitting us to identify them. 46 Moreover, we can explain the 
spin-polarization as due to the different spatial extend of the molecular orbitals in each of the 
spin channels and the corresponding overlap with the electrodes. Hence, the spin-polarization is 
rather an effect of the geometry of the molecular orbitals at play rather than due to a spin-polarized 
density of states. 

Projecting the density of states onto atomic orbitals is also instructive. Figure [3] (c) depicts 
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the PDOS onto the atomic TM-d and Carbon-p states. This permits us to corroborate the above 
conclusion. Indeed, we can see that while the HOMOf has a large component on Carbon-p states, 
the HOMO), is basically a TM-d orbital. This same conclusion, but for different orbitals, is deduced 
from the LUMO composition. 47 We can then conclude that the larger contribution to the electronic 
current of the majority spin (j') orbitals is due to the contribution of Carbon-p states, and hence of 
the ^-orbitals of the Cp ligands revealed in Fig. [2] 

Transport in the contact regime 

To mimic the contact regime, we approach the right electrode to the molecule at a distance of 
6?2 =2.57 K(d\ —2.12 A ). 48 The present electrode configuration has been chosen so as not to exert 
any pressure on the molecule. Hence, the electrodes induce negligible distortions of the molecular 
geometrical parameters. The charge transfer to the molecule is still small (0.177 e) and magnetic 
moment of the molecule reaches 1.895 fig (with 0.478 fls, 1.630 / 1 b, -0.234 /t# and 0.021 /is for 
Fe, Co, C and H, respectively, which are similar to the ones described in Table 1). Overall, this 
analysis and the one described for tunneling conditions indicate a modest effect of the electrodes 
and vdW-forces on the properties of the molecule for the two-probe system. 

Figure [4] (a) shows the transmission at zero bias as a function of the electron energy. As 
expected the transmission is larger than the transmission in the tunneling regime, leading to a 
total conductance at the Fermi level of G(£>) = 0.073Go. At the Fermi level, the majority spin 
channel exhibits a transmission probability one order of magnitude higher than the minority spin 
channel. Thus, the molecule maintains its spin-filter character (SFE — 86%). 

For both spin channels, eigenchannel analysis— shows that two scattering states provide the 
major contribution to the transmission function in the whole energy range. In particular, the two 
most contributing scattering states provide very similar contributions at the energies corresponding 
to peaks PI, P2, P3 and P4 (see Fig.|4jb,d)). Although, it is difficult to identify eigenchannels by 
visualizing them, 50 the perfect energy alignment between these peaks and the ones observed in the 
PDOS (Fig. |4](c)) allows again to assign PI (P2) and P3 (P4) peaks to transmission trough HOMOj^ 
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(HOMO|) and LUMOf (LUMO|) molecular orbitals. Moreover, assuming Breit-Wigner-like res¬ 
onances for the transmitting MO, we have fitted the corresponding transmissions with Lorentzian 
functions 41 which permits us to confirm that the peaks in the transmissions nicely corresponds with 
the molecular levels in the PDOS (see Supporting Information). Moreover, the Lorentzian fitting, 
albeit imperfect, shows that the LUMO transmission dominates at the Fermi energy for both spin 
channels. As the electrode approaches the molecule, the contributions of the LUMOs grow, with 
no reversal of molecular character in the electron transmission. 

As can be seen, the 4 frontier-orbital peaks shift to lower energies with respect to their en¬ 
ergy position in the tunneling regime. This is due to the enhancement of the molecule-electrode 
interactions which also induce a more pronounced broadening of the involved molecular levels. 
Interestingly, the largest hybridization is observed for PI where the broadening increases roughly 
a factor of 6. This behavior can be traced back to the larger overlap of the HOMO^ with the 
approaching electrode. 

Our transport calculations carried out for the staggered conformer show that the symmetry of 
the molecule does not affect the spin-filter character of the molecular junction (see Supplementary 
Information for more details.) 

Finite-bias results 

The above results imply that FeCoCp 3 is a good spin filter in the linear-response regime. In this 
section, we go beyond the linear-response regime. We computed the electron current for both spin 
channels (I|,I|) as a function of the applied bias using eq. ??. 

Figure|5£a) shows the electron transmission that enters the Landauer equation, eq. ??, evaluated 
for three different bias. The upper panel shows the majority spin transmission. We find that the 
HOMO-LUMO gap increases with bias and the transmission in between the two peaks decreases. 
For the minority spin (lower panel) the bias effect is negligible. Overall, the effect of the bias is 
small and using the zero-bias transmissions seems justified. Nevertheless, it is interesting to both 
understand why the bias effect is small and why the effect is not noticeable for the minority spin. 
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The effect is small because the molecule is basically bound by dispersion forces, hence the 
molecular electronic structure presents small perturbations from the electrodes. The presence of 
an external electrical field acts on the polarization of the molecule. Here, the fields are so small that 
this effect is negligible. The flowing of a current through the molecule is a larger effect, leading to 
a change in the steady-state charge of the molecule. However, the HOMO stabilizes by trapping a 
very small amount of charge and in the same degree the LUMO empties, contributing to an almost 
zero change in charge state. This leads to a small opening of the HOMO-LUMO gap. 

As we have previously seen, the minority-spin molecular orbitals are less coupled to the sub¬ 
strate. Hence, the opening of the HOMO-LUMO gaps is negligible. Interestingly, the effect of the 
bias on the magnetic moment of the molecule is negligible (i.e it goes from 1.895 /is to 1.860 /is 
when the bias increases 0.5 V, see Tables 3-4 in the Supporting Information). 

Figure [5fb) shows the current computed using the Landauer equation, eq. ??. In the linear 
regime the majority-spin current, 1^ is appreciably larger than I| due to the higher conductivity for 
majority than for minority spin channels (0.0675Go vs 0.0055Go). Such a large difference between 
I| and I| is still observed as the voltage further increases. As a result, FeCoCp 3 acts as a spin filter 
in the whole bias voltage range with a large current polarization, CP ~ 84 %. The inclusion of bias 
in our calculations does not change the conclusion that this molecule is an excellent spin filter with 
a CP close to 90%. 

These results are in contrast with the ones of Ref. 51 where they find that the CP in a Fe-C 7 oC 7 o- 
Fe junction goes from 78% at zero bias to 20% at 0.5V. The large difference between our results 
and theirs can be traced back to the very different interaction of the molecules with the electrodes. 
While in our case the molecule is physisorbed by van der Waal forces, in their case, a strong 
covalent interaction rules the charge flow through their Fe-C 7 oC 7 o-Fe junction. 

Spin flip effects 

The electronic current can yield energy to the molecular spin degrees of freedom, and hence change 
the spin state of the molecule. As a consequence spin excitations can reduce the spin-filtering 
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capabilities of the device if the excited states corresponds to different spin alignments. Let us 
briefly describe spin transport through the FeCoCp 3 molecule. 

The MAE of the molecule is 1.64 meV as described above. The molecular axis aligning the Fe 
and Co atoms is a hard axis. Hence, the molecular ground state corresponds to a spin of 1 in the 
easy plane described by the Cp ligands which corresponds to a S z — 0 if the molecular axis is taken 
as the z-axis. In this conditions, the electron spin is contained in the molecular easy plane. As we 
have seen before, the spin-polarization with respect to an axis on this plane will be very large, well 
above 80% in all the cases analyzed above. Precession of the spin-polarization axis will be small, 
and the spin current will be polarized in an arbitrary axis contained in the molecular plane. 

A spin Hamiltonian can be written that reproduces the MAE for this 5=1 molecule. We can 
easily see that 

H spin =DSl (7) 

In the present case the value of D is 1.64 meV. From here we see that the first excitation is indeed 
equal to 1.64 meV and it corresponds to flipping the spin from 5, = 0 to \S Z \ — 1, i.e. from the easy 
plane to the hard axis. Hence, electrons with energy above the first-excitation threshold (biases 
above 1.64 mV) can flip the molecular spin out of the easy plane if they flip their spin. A simple 
calculation 19 shows that the incoming electron has a probability of 1/3 to flip its spin in the present 
case. As a consequence the CP goes from a value close to 100% to 33% when the absolute value 
of the applied bias goes above 1.64 mV (in the case where the intrinsic spin-polarization due to 
the electronic structure is the 84% of the previous section, the spin polarization above the spin-flip 
threshold becomes 28%). 

This description is valid both for the tunneling and the contact transport regimes, since only 
the molecular MAE and spin multiplicities enter it. 

Thermoelectric effects 

Motivated by the different ratios x a /x a at the Fermi level for minority and majority spin channels, 
we evaluate whether a spin-polarized thermopower current can reduce the spin polarization of the 
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total current. This is of importance because the spin filtering capacities may not be maintained 
in the presence of a temperature drop (AT — T L — T R ) across the junction. This physical situation 
can be reached when the electrodes are contacted in a different way, and current dissipation in the 
electrodes may lead to different temperatures. 

For this purpose, we take a temperature drop AT — — 1 OK between electrodes and compute 
the spin-polarized electron current, I a with o =j\J, using eq. ??, for different Bias. The current 
polarization, CP (eq. ??), obtained in each case as a function of the average electrode temperature 
T (T — ( Tl + 7)0/2) is shown in Fig. [(Ja). For the sake of comparison, we plot the CP values 
obtained when both electrodes are at exactly the same temperature. 

At zero (Bias = 0) and extremely low bias (Bias — 2 x 10 6 V), we see that thermal effects 
induce a drop of the CP value from ~86% when both electrodes are at the same temperature to 40- 
50% in presence of a small temperature gradient. However, the excellent spin-filtering capabilities 
are restored as soon as the bias voltage is slightly increase; the CP reaches again 86% when the 
bias is 0.02 V. 

To understand such thermal effects on the current polarization, one simply needs to make use 
of the linear-response limit of the spin polarized electron current which tells us that l a = I x a pl'jj — 
G a V + G a S a AT ; a =j\ j, (see eq. ??). From this expression, we can clearly establish two limiting 
behaviors: one dominated by thermoelectric current 7^' at low biases, and the other one by the bias, 
I x 0 . when the bias becomes larger than a critical bias, V c , given by V c ~ k R AT. 

Let us focus on the first case, Fig.[6tb), where I a can be approximated by G a S a AT (a =1\D- 
Here, CP is reduced to (G^S^ — G^Sj.) / (G^Sj + G|S|). Hence, the Seebeck coefficient (Sa) times 
the conductance (G a ) for the two spin channels are the key ingredients of the current polarization. 
The spin-dependent conductances G| and G; with average values 5250 nA/V and 409 nA/V, 
respectively, barely change in the studied temperature window. In addition, the spin-dependent 
Seebeck coefficient as a function of the electrode temperature plotted in Fig. [6fc) shows that \S^\ 
is roughly four times lower than |SJ. Asa result, |G^5|| is crudely three times larger than | G|SJ 
(see Fig. Hid)) which explains the 40-50 % of current spin polarization observed in Fig. [6la). 
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With regards to the second case where I G ~ G a V, the current spin polarization is here simplified 
to CP ~ (Gj — G|)/(Gj + G|). Therefore, the excellent spin-filtering capacities (CP = 86%) found 
in this case can be traced back to a much higher conductance for majority than minority spin 
channels. 

Summarizing, thermoelectric effects in this type of molecular junctions lead to a strong sup¬ 
pression of the otherwise excellent spin-filtering properties of the molecules when the electronic 
transport is governed by the thermoelectric current. 

A different thermal effect is the one given by a homogeneous temperature. As the temperature 
rises, the direction of the molecular spin can change. Indeed, at ~20 K, the ambient temperature 
is large enough to induce spin flips, similar to the spin-flips we have described in the previous 
section. 


Summary and Conclusions 

Using DFT calculations together with a NEGF implementation of electronic transport equations, 
we have evaluated the gas-phase, adsorption and transport properties of a CoFeCp 3 . The motivation 
to do so is the spin (S=l) of the gas-phase molecule, and its magnetic anisotropy (MAE=1.64 meV). 
These two properties are good characteristics for a tentative molecular-based spin filter. 

The molecular spin is largely localized on the Co atom, and the Fe atom is basically not mag¬ 
netic. This is due to the charge transfer originating in the Cp ligands, and is in agreement with 
what is found for cobaltocene and ferrocene. 

On a Cu(l 11) surface, we find that the molecule binds via dispersion forces and that the charge 
transfer is negligible, hence keeping the above molecular properties. The molecules present two 
conformers, one where the Cp rings are aligned, eclipsed conformer, and a second conformer where 
the Cp are alternatively rotate in a staggered fashion. We find that systematically the eclipsed 
conformer is more stable. 

The transport properties of the molecules are computed in the tunneling and contact regimes. 
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On the adsorbed-molecule setup, a second electrode is approached. We have used an electrode¬ 
molecule distance of 5.15 A to characterize the tunneling regime. The contact regime corresponds 
to a molecule-electrode distance of 2.72 A. We find that the Fermi energy is in the middle of the 
HOMO-LUMO gap and that the transmission is largely dominated by the tail of the LUMO res¬ 
onance. Due to the large contribution of the Cp ligands to the majority-spin HOMO and LUMO 
we find a large electron transmission for the majority spin channel. At the same time, the electron 
transmission through the minority-spin channel is smaller due to the prevalence of the TM-d or¬ 
bitals. As a result, we find a strong spin polarization in the current, with a polarization of 98% in 
the tunneling geometry and 86% in contact. 

When voltage is applied across the molecular junction, we find a small opening of the HOMO- 
LUMO gap in the majority-spin channel, while a negligible effect for the minority-spin one. The 
current spin polarization is very constant, changing from the above 86% at 0 V to 83% at 0.5 V. 
The behavior with bias is very weak due to the weak coupling of the molecule to the electrodes and 
the negligible charge transfer. However, as the bias increases inelastic channels open that further 
reduce the spin polarization of the current. 

For biases larger than 1.64 mV, equivalent to the MAE of the molecule, electrons can flip the 
molecular magnetic moment out of the easy plane. As a result the spin of electrons also change 
and the spin polarization is reduced. For the first excitation threshold this reduces the current 
polarization to 33%. 

Also thermoelectric effects in the absence of applied bias lead to a strong suppression of the 
otherwise excellent spin-filtering properties of the molecules. When bias is applied, the much 
larger bias contribution overrides the small thermopower and the spin-filtering properties of the 
molecular junction are recovered. 

In conclusion, a superficial analysis of our calculations would show the triple-decker molecule 
CoFeCp 3 as an excellent current spin filter. However, spin-flip processes and thermocurrents have 
very negative consequences for this type of device. A negligible temperature difference between 
electrodes can rapidly diminish the spin-filter efficiency when the electronic transport is governed 
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by thermoelectric currents. Moreover, ubiquitous spin-flip inelastic effects need to be considered 
when evaluating the spin-filtering properties of a molecular junction. 
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Figure 2: Doubly-degenerated frontier molecular orbitals for the D^h conformer. Plotted isovalues 
are 10 % of the maximum ones. Red (green) indicates positive (negative) values of the real part of 
the wavefunction. 
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Figure 3: (a) Electron transmissions from the left to the right electrode as a function of electron 
energy referred to the Fermi energy. (b)Computed spin-polarized PDOS onto frontier molecular 
orbitals, (c) Computed spin-polarized PDOS projected onto TM-d and Carbon-p atomic orbitals. 
Tunneling regime d2=5.15 A. Total magnetization^ 2 jig. Mulliken charge transfer to the 
molecule = 0.098 e 
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Figure 4: Contact regime cb ~ 2.57A. (a) Electron transmissions from the left to the right 
electrode as a function of electron energy referred to the Fermi energy, (b) the PDOS onto TM-d 
and Carbon-p states, (c-d) Eigenchannels contribution to the average transmission for majority 
and minority spin channels (the color-code only denotes the ordering of the transmission 
eigenvalue at each electron energy and may not be connected to the symmetry or other common 
properties of the eigenchannels).Total magnetization^ 1.92/1/;. Mulliken charge transfer to the 
molecule = 0.177 e. 
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Figure 5: (a) Electron transmission as a function of electron energy for 0, 0.3 and 0.5 V applied 
between the two metallic electrodes. The majority spin (upper panel) shows that the 
HOMO-LUMO gap is broadened as the bias is increased and there is a small drop of transmission 
at the Fermi energy. For the minority spin (lower panel) the shifts are negligible and there is a 
very minor bias dependence, (b) I-V characteristics in nA and V. The SFE is shown in brackets. 
Dot lines indicates the results obtained with the linear approach l a — G a V 
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Figure 6: (a) Current polarization (CP) as a function of the average electrode temperature T (i.e. 
T = (Tl + Tr)/2) for different Bias when AT = — 10/C The black dash line indicates the CP 
values obtained when both electrode are at exactly the same temperature, (b) Majority (black 
lines) and minority (red lines) spin polarized currents for an applied bias of 2 x 10 6 V 
(AT — —\0K). Dash and dash-dot lines indicates 1^ and l T a h contributions to the spin-polarized 
current I a shown in solid lines; o =|. | (see text), (c) Spin dependent Seebeck coefficient S a 
(/iV/K) as a function of the temperature (K). (d) The Seebeck coefficient times the linear 
conductance (nA/K) for the two spins as a function of the temperature (K). 
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